hep-ph/OlU^UgB 
CERN-TH/2001-200 
IFP-796-UNC 
MPI-PhT/2001-27 



Measurement of the Lifetime Difference of Mesons: 

Possible and Worthwhile? 



A. S. Dighe^'f], T. Hurth^f], C. S. Kim ^f|, and T. Yoshikawa 

a: Max-Planck-Institut fiir Physik, Fdhringer Ring 6, D-80805 Miinchen, Germany 
b: Theory Division, CERN, CH-1211 Geneva 23, Switzerland 
c: Department of Physics and IPAP, Yonsei University, Seoul 120-749, Korea 
d: Department of Physics and Astronomy, University of North Carolina, 
Chapel Hill, NC 27599-3255, USA 

(February 1, 2008) 
Abstract 

We estimate the decay width difference ^Td/Td in the Bd system including l/m^ 
contributions and part of the next-to-leading order QCD corrections, and find it to 
be around 0.3%. We explicitly show that the time measurements of an untagged 
decaying to a single final state isotropically can only be sensitive to quadratic terms in 
AFrf/r^, and hence the use of at least two different final states is desired. We discuss 
such pairs of candidate decay channels for the final states and explore the feasibility of 
a AF^/r^ measurement through them. The measurement of AF^^ would be essential 
for an accurate measurement of sin(2/3) at the LHC. The nonzero width difference can 
also be used to identify new physics effects and to resolve a twofold discrete ambiguity 
in the B^-Bd mixing phase. We also derive an upper bound on the value of AF^/F^ 
in the presence of new physics, and point out some differences in the phenomenology 
of width differences in the Bs and Bd systems. 
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1 Introduction 



Within the standard model (SM), the difference in the decay widths of Bd mesons is CKM- 
suppressed with respect to that in the Bg system. A rough estimate leads to 

^~^-A=«0.5%, (1) 

i- d i- s 

where A = 0.225 is the sine of Cabibbo angle, and we have taken AFs/F^ ~ 15% |1|, ^ |^. 
Here Td(s) = (Xl + ^h)/'^ is the average decay width of the light and heavy Bd(s) mesons 
{Bl and Bh respectively). We denote these decay widths by F^, Th respectively, and define 
AFrf(s) = Fj^ — Th- No experimental measurement of AF^ is currently available. Moreover, 
no motivation for its measurement (other than just measuring another number to check 
against the SM prediction) has been discussed, and hence the study of the lifetime difference 
between B^, mesons has hitherto been neglected as compared to that in the Bg system. The 
phenomenology of the lifetime difference between Bg mesons has been explored in detail in 

Si- 

With the possibility of experiments with high time resolution and high statistics, it is 
worthwhile to have a look at this quantity and make a realistic estimate of the possibility 
of its measurement. At LHCb for example, the proper time resolution is expected to be as 
good as At ^ 0.03 ps. This indeed is a very small fraction of the B^ lifetime {tb^ ~ 1.5 
ps [0), so the time resolution is not a limiting factor in the accuracy of the measurement, 
the statistical error plays the dominant role. Taking into account the estimated number of 
Bd produced — for example the number of reconstructed B^ ^ J/ip Ks events at the LHC 
is expected to be 5 x 10^ ([0 table 3) — the measurement of the lifetime difference does 
not look too hard at first glance. Naively, one may infer that if the number of relevant 
events with the proper time of decay measured with the precision Ar is A^, then the value of 
AFd/Frf is measured with an accuracy of 1/^/N. With a sufficiently large number of events 
A^, it should be possible to reach the accuracy of 0.5% or better. 

The measurement of ATd/Td is in reality harder than what the above naive expectation 
may suggest, since most of the quantities that involve the lifetime difference are quadratic 
in the small quantity ATd/Td- In fact, as we shall explicitly show in this paper, the time 
measurements in the decays of an untagged Bd to a single final state are sensitive only to 
(AFrf/Fd)^. This implies that in order to discern two different lifetimes, the measurements 
need to have an accuracy of (AF^/Fd)^ ~ 2.5 x 10^^, which is beyond the reach of the 
currently planned experiments. 

However, the combination of lifetimes measured in two different untagged decay chan- 
nels may be sensitive to linear terms in AF^/F^. We explore three pairs of such untagged 
measurements in this paper: (i) lifetime measurements through decays to self-tagging (e.g. 
semileptonic) final states and to CP eigenstates, (ii) CP even and odd components in the 
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decay mode — > J /ipK*{KsTi^), and (iii) time-dependent untagged asymmetry between 
Bd ^ J/^Ks and 5, ^ J/^Kl. 

The conventional "gold-plated" decays for f3 measurement, J/ipKs and JjipKL, neglect 
the lifetime difference while determining sin(2/3). For an accurate determination of (3, the 
systematic errors due to AF^/Frf need to be taken into account. Moreover, there is the 
possibility that the measurement of the lifetime difference leads to a clear signal for new 
physics. Furthermore, if the lifetime difference is neglected, the ambiguity (3 ^ (7r/2 — /5) 
remains unresolved. Observables that are sensitive to the lifetime difference may resolve this 
discrete ambiguity under certain conditions. 

The observables mentioned above can also give an independent measurement of cos(2/3) 
in principle. In order to be able to do this, however, the theoretical uncertainties on AF^ 
need to be minimized. Therefore, we start by presenting in Sec. ^ a detailed calculation 
of AFrf, including 1/mf, contributions and part of the next-to- leading order (NLO) QCD 
corrections. The NLO precision in the width difference AF^ is also essential for obtaining 
a proper matching of the Wilson coefficients to the matrix elements of local operators from 
the lattice gauge theory. 

The rest of the paper is organized as follows. In Sec. ^ we explicitly demonstrate the 
quadratic dependence on AF^/F^ of quantities measurable through untagged B decays to a 
single final state. We explore the combinations of decay modes that can measure quantities 
linear in AF^/F^. We calculate the corrections due to AF^ to the measurement of sin(2/3) 
through Bd —>■ J/ipKs, and also indicate the possibility of the AF^ measurement through 
tagged decays to CP eigenstates. In Sec. ^, we point out important differences in the upper 
bounds on AF^ and AF^ in the presence of new physics, and elaborate on the possibility of 
detecting new physics and resolving discrete ambiguities in the mixing phases through them. 
We summarize our findings in Sec. ^. 



2 Estimation of AF^^ 
2.1 Basic definitions 

We briefly recall the basic definitions: in the Wigner-Weisskopf approximation the oscillation 
and the decay of a general linear combination of the neutral flavour eigenstates B^ and Bd, 
a\Bd) + b\Bd), is described by the time-dependent Schrodinger equation 

'10 -(--HO- 

Here M and T are 2x2 Hermitean matrices. CPT invariance leads to the conditions 
Mil = M22 and Fn = F22. Exact CP invariance would imply M21 = M12 and F21 = F12 
(a phase choice, namely CV\Bd) = ~\Bd) ,CV\Bd) = —\Bd) is made). Independent of the 
choice of the unphysical phases, CP invariance (in mixing) would imply Im(M2iF2i) = 0. 
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The mass eigenstates, the hght Bl and the heavy Bh, are given by 



\BL,H)=p\Bd)±q\B,) (3) 

with the normahzation condition \q\^ + |pp = 1. Only the magnitude \q/p\ is measurable, 
the phase of this quantity is unphysical and can be fixed arbitrarily by convention. 

The mass difference and the width difference between the physical states are defined by 

Am^Mn- Ml, AT = Tl-Th , (4) 

such that Am > 0, AF^ > in the SM. The real and imaginary parts of the eigenvalue 
equations are the following: 

(Amf-^(Ar)2 = (4|M2i|2-|r2i|2), (5) 
AmAr = -4Re(M*ir2i). (6) 

With the help of the CP- violating parameter S 

-2Im(M|,r,i) 

The effect of CP violation due to mixing on the mass difference Am and on the lifetime 
difference AP may be explicitly shown: 



,2 _ 4|M2i_p-f |r2i'2 

1 + 52 



(Am)^ = ' , (8) 



In the limit of exact CP invariance (d = 0) the mass eigenstates coincide with the CP 
eigenstates, CV\Bh) = —\Bh) and CV\Bi) = +|-Bl) and the mass difference and width 
difference are given by Am — 2|M2i|,Ar = 2|r2i|. However, even with a non-zero S, 
taking into account that 5 is constrained by the upper bound \5\ < |r2i|/(2|M2i|) and 
r2i/M2i = 0{ml/mj), we can write 



Am = 2IM21I 



.4' 



(10) 



IM21I 

We shall neglect the terms of 0{ml/mf) ~ 10~^ in our calculations. 
2.2 Method of calculation 

In the following we consider the two off-diagonal elements M21 and r2i, which correspond 
respectively to the dispersive and the absorptive part of the transition amplitude from B^ 
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Figure 1: Schematic representation of Feynman diagrams for M12 = and 
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to Bd- We follow the method of 



which was used there in the Bg-Bg system (see also 



Within the SM the well-known box diagram is the starting point of the calculations. M21 
is related to the real part of this diagram (see Fig. 0). The important QCD corrections are 
most easily implemented with the help of the standard operator product expansion. Because 
of the dominance of the top quark contribution, M21 can be described by a local Ai? = 2 
Hamiltonian below the mw scale: 
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{Bd\H^Ff "^{Bd) 



O 



m 



w , 



n 



AB=2 
eff 

Q 



G 



^ iV*VtdYC'^{mt,mw,f^)Qi^) + ii.c. , 



l_67r2 

(bidi)v 



(11) 

(12) 
(13) 



The Wilson coefficient C"^ contains the short- distance physics. It is known up to NLO 
precision []TU[. The hadronic matrix element {Bd\Q{fi ~ 'mb)\Bd) will be discussed below. 

In the standard model, r2i is related to the imaginary part of the box diagram. Via the 
optical theorem it is fixed by the real intermediate states. Therefore, only the box diagrams 
with internal c and u quarks contribute (see Fig. [I|). In contrast to the Bg-Bg case where 
the intermediate cc contribution is the dominating one, because of its CKM factor (V],^Ks)^, 
over the uu, the cu and the uc contribution (see Sec. |0| ), in the Bd-Bd all four contributions 
have to be taken into account. In the effective theory where we integrate out the W boson, 
r2i is given by: 



1 
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2Mb, 



{Bd\lm I d^xTH 



AB=1. 
eff ' 



X] 



^eff 



Bd) 



(14) 



where 



AB=1 
eff 



^(KT^Kd E C.Qr + v:,Vud E C^QT + v:,v,d e c^qt + 



=1,2 



=1,2 



=1,2 



+ v^v^d E c^QT - v:,Vu E cm- 

i=l,2 i=3 



(15) 
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The operators are denote color indices) 

Qf = (Piqj)v-A{q'jdi)v-A, 

Qs = {hdi)v-A{(ij(lj)v~A-, 

Qb = (bidi)v-Aiqjqj)v+A, 
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Q4 



{biqi)v-A{q'jdj)v-A, 

{bidj)v-A{qjqi)v-A, 
{hidj)v-A{qjqi)v+A- 



(16) 
(17) 
(18) 



The penguin operators - Q& have small Wilson coefficients and are therefore suppressed 
with respect to the four-quark operators - which all have the same two Wilson coefficients 
Ci and C2. In the leading logarithmic approximation we have: 



'a{Mw)' 
. a(/i) . 



2/3o 



W) 



(19) 



where /3o = (llA^-2/)/3 = 23/3 and = ±6(1±A^)/A^. The coefficients to NLO precision 
can be found in ||lT| . 

Because there is another short-distance scale, the bottom quark mass, the operator prod- 
uct of two AB = 1 operators can be expanded in inverse powers of the bottom quark mass 
scale in terms of local AB = 2 operators: 
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B, 



f/r(o)|5.) 



„ ml 



(20) 
(21) 



These matching equations fix the values of the AB = 2 Wilson coefficients En- The corre- 
sponding four quark operators 0„ are the following: The operators Q and Qs, 



Q — {hdi)v-A{i>jdj)v-A-, 
Qs = (bidi)s-p(bjdj)s-p, 



(22) 
(23) 



represent the leading order contributions. Their matrix elements are given in terms of the 
bag parameters, B and Bs, the mass of the Bd meson Mb^, and its decay constant /b/- 



{Bd\Q\Bd) 
[Bd\Qs\Bd) 



/^M|2^ + ^ 



B, 



-flMl 



Ml 

Bd 



2N 



N 



-B. 



(24) 
(25) 



In the naive factorization approximation, B and Bs are fixed hj B = Bs = 1. Reliable 



lattice calculations for B and Bs are already available [0. We note that to NLO precision 
one has to distinguish between the pole mass and the running quantity 

,2 



mf,(yu) = rrih 



1 - 



a. 



71 



. mi 3 , 



(26) 
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using the MS scheme. 

The l/mfe corrections are given by the operators 



Ri 
R2 
R3 

Rq 



— {bidi)s-p{bjdj)s+p, 

1 



mt 



mt 



(6,(l-75)zD^rf.)(6,Y(l-75)f^i) 



Qs + ^Q + Qs, 



(27) 
(28) 
(29) 
(30) 
(31) 



where Qs has the "interchanged" color structure as compared to Q^. There are also "color- 
interchanged" operators Ri and Q corresponding to Ri and Q. We note that these I /nib 
operators are not independent, the relations between them are in fact the equations of 
motion. 

The matrix elements of these operators within the Bg-Bs system were estimated in [0 
using naive factorization, which means that all the corresponding bag factors were set to 1. 
For the Bd-Ba system the analogous results are: 



{Bd\Ro 
{Bd\Ri 
{Bd\Ri 

{Bd\R2 

{Bd\Rz 
{Bd\R^ 
{Bd\RA 
{Bd\Ri 



\Bd) 
\Bd) 
\Bd) 
\Bd) 
\Bd) 
\Bd) 
\Bd) 
\Bd) 



\ N 



2 ,,2 rrid 2N + 1 
^''rub N 



f B^ ^B^ 
flXB. 
flA 



'Ml 



mi 
. ml 

. ml 
(Ml 



^ - 1 



- 1 



ml 



= 0, 

, 

1 - 

2N + 1 

2N 
N + 2 

2N ' 



■{Bd\R2\Bd) 



mt 



-SI A {% - \ 



1 



(32) 

(33) 
(34) 

(35) 
(36) 
(37) 
(38) 
(39) 



We neglect terms proportional to md/mb, the other terms proportional to (M'^^/ml) — 1 are 
of order KqcD/mb- 

In the matrix elements {Ri) (eqs. (P^-(P^)), we use the pole mass mb- There is a 
subtlety involved here: as discussed in there are terms of order as and of leading power 
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in rrib in the matrix element of Ro to NLO precision. In view of the relation (PT]), it is not 
surprising that there are such terms. In the scheme - which was used in ^ and which is also 
used here - these terms are subtracted in the matrix element (Ro) while taking into account 
the leading NLO contribution. Then the (Rq) matrix element is still of a subleading nature. 
The specific subtraction scheme for the factorized matrix elements (Ri) corresponds to the 
use of the pole mass in eqs. (|5^-(|5UD. Of course this specific choice for the matrix elements 
has to be taken into account if the NLO results are combined with a lattice calculation of 
the {Ri). 

There is an additional remark in order. We estimate by the cut of the partonic 
diagrams. The underlying assumption of local quark-hadron duality can be verified in the 
Bg-Bg system, in the simultaneous limit of large N and of small velocity H, therefore one 
expects no large duality violations. In the B^-B^ system the small velocity argument fails 
since the uu, uc and cu intermediate states contribute significantly, and the larger number 
of light intermediate states leads to a larger energy release. We follow ref. |^ and make the 
assumption that the duality violations in the Bd — Bd system are also not larger than 10%. 
In order to test this assumption one should include all corrections up to that accuracy. 

2.3 Analytical results 

In this section, we present an analytic expression for including 1/mf,, penguin and part of 
the NLO corrections. If one takes into account the error inherent in the naive factorization 
approach to the matrix elements of the subleading operators R, it seems to be a reasonable 
approximation to keep at least all terms up to an accuracy of 10~^ -p^eading^ keep also 
higher order terms in order to check the accuracy of our approximation. 

In the effective theory of the AB = 2 transitions the matrix elements of the l/rrib op- 
erators (R) are formally suppressed by a factor of the order of 0.1 with respect to those of 
the leading operators Q and Qs- The natural variable z = m'^/ml also formally introduces 
a suppression factor of approximately 0.1. The NLO contribution has formally an extra 
suppression factor {ag/ATr) of order 0.01. Within the effective theory of the AB = 1 Hamil- 
tonian, the combination K' = (^peng^-^dom ^^^^ ^// _ (jpeng^peng g^j^g suppressed by almost a 



factor 0.01 and 10~^ respectively, with respect to the combination K = C'^°"^C'^°'^, where 
(jpeng (jg^otes the Wilson coefficients of the penguin operators Qs-.-Qe and C**^"™ that of the 
dominating operators Qi'(2)- The contribution due to K" therefore can be safely neglected. 
Schematically our analytical result for has the following form: 

r2i = K{Q) (40) 

+ K{R)i0il) + 0iz) + 0iz^) + {0iz^)}) (41) 

+ K'{Q)iOil) + 0{z) + {Oiz')}) (42) 

+ K' {R) {0(1) + {0{z')} (43) 
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+ «,/(47r) i^(Q)0(l) 



(44) 



where Q represents the leading order operators Q and Qs- The terms inside the curly 
brackets are the ones that we calculate only to estimate the errors. In the presentation of 
the results the following combinations of the Wilson coefficients are used: 



K' 



2(3CiC3 + C1C4 + C2C3), K'^ = 2C2C4, 
2(3CiC5 + CiCq + C2C5 + C2Cq), 



C1C2 



(45) 
(46) 
(47) 



and the common factor of [—G'j;'ml/(24:7TMBj)] is implicit in the following equations (^8[) , 

(1), (0), (H). 

In the leading log approximation we calculate the Qf and the contributions to r2i. 
By extracting the absorptive parts of the cc, uc, cu and uu intermediate states, we can find 
the off-diagonal element. For this leading contribution (|^), after replacing V*^Vud by the 
unitarity relation, we get to all orders in z: 



r leading 
21 



K, + -K2){Q) + {K,-K2){Qs) 



+ iV:M{V,lVtd) [{3z{Ki + K2) - 3z^K2 - z\Kr - K2)){Q) 
+{6z\K^ - K2) - ^z\K, - K2)){Qs) 



Ki + -K2 



iQ) 



z{K^ + 2^2 
+71^(1 + 2z){K^ - K2){Qs) 
+ f (1 - zf [-{2K^ + K2) + z{K2 - K,)] + 



Ki + -K2 



(Q) 



+ 



;i - z)^{2 + Az){K2 - Ki) - {K2 - Ki)] {Qs)} 



(48) 



This result, (^Bj), confirms the findings of 0. 

In the equations (P^),(pUD, (pT|), and ( |52D we give the penguin and l/rrib contributions. 
The coefficients of (VJ^^td)^ in these results can be checked with the z — > limit of the results 
in the literature for the Bs syatem. In this sense our results are consistent with the findings 
in the Bs system given in 0. 

The l/mb corrections to the operators Qf and Q'2' give [see the term (^)] 



pl/rUb 
J- 21 



-2[K,- -K2 ] {R2) - 2K,{R,) + 2K2{R,) 



+ {V:,V,d){V,lVu)[-l2z'K,{{R,) -2{R,)) 

+ 6z^K2{{R2)+MR3) + 2{R,)) 

+ {k{R)0{z% 



(49) 
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The term in curly brackets in (H^) can be written as 



{...} = {v:^v,,){v:,vu) ^ 

[12/iri(2(i?i) - (i?2) - 6(i?3)) - 12z3ir2((i?2) + 6(i?3) + 2(i?4))] 



(50) 



The penguin contributions [terms 



(O] are 



+ (-2(i?2) - 2{R^))K[ + ((i?2) + 2(/?4))i^2 

+ (v;;v;rf)(v;*,\4d)(3zir( + szir^ - ?>zK'^{q) 

+{K'{Q)0{z^)+ K'{R)0{z)} , 



where the terms in curly brackets (and the lower order ones) may be written as 



{...} = {v:,v^){v:,vu) ^ 

-{l2z^R,)-2Az^R,))K[ 

+ {Qz^R2) + 24z^{R3) + 12z^R^))K'^ + Az^R2)K'^ 



(51) 



(52) 



The NLO QCD correction T 
the limit 2 — > of their results^: 



NLO 
21 



«s/(4vr) K (Q) [term (Hj)] is found from @] by taking 



^NLO 

21 



109 



248 



V 157^ 
3 ^ 18 , 



iQ) 



+ 



(lOK, + '^K, + Ik,) In (^) - (34^4 + 10K2)ln 
\ 3 3/ \mhJ Vmft 



+ 



40^^ 248^^ S-n 
3^9 ^ V 3 

327^3 - y + y i^2 ) In 



rrib 
128\ 

/i2 



(Q) 



27 9 Knib 



+ 



ICs 
3 Co 



- (167^4 + 167^2) In 
mQ)-8{Qs))} . 



{Qs) 



(53) 



The explicit //i and fi, dependence in ( |5B| ) cancels against the /i dependence of the Wilson 
coefficients of the hamiltonian Ti^^"^ (|T3p and the fi dependence of the matrix elements of 
the AB = 2 operators at the order in as we take into account. For a proper matching with 



^ We add only the leading contribution of the NLO QCD corrections for the term (Vj^Vfd)^- The lead- 
ing terms of the contributions for the terms {V*jjVcd){Vi*^Vtd) and {V*jjVcd)'^ cancel out through the GIM 
mechanism. 
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lattice evaluations of these matrix elements it is important to note that the results in 
are based on the NDR scheme, with the choice of 75 and the evanescent operators as given 
in eqs. (13)-(15) of ||. 
The net is 

■p Y^leading . -pl/m,;, j^peng . -r^NLO /'C/tN 

i 21 — i 21 + J- 21 + J- 21 + J- 21 ) l^4j 

with the implicit multiplicative factor of [— G'|'m^/(247rMB^)]. 
2.4 Numerical results 



Let us now calculate the numerical value of AF^;. From eq. ([lOD , AF^ can be approximately 
written as 

AFrf ^ -2|M2i|Re^ = -Am Re^ . (55) 
' ' M21 M21 ^ ^ 

where M21 [see eq. (0)] is given by 

= 7TWS^^(^tlVuySo{x,){Q). (56) 
(47r)^(2MijJ 

Here = rh'f/M^, rjB is the QCD correction factor and 5*0 is the Inami-Lim function: 

Using the results obtained in the previous section, we can write down the width difference 
(normalized to the average width) in the form 



AF\ /Am 



V F Jb, V F /, 



/C X 



rtt I ^ ^sMl 3 1 



(5J 



The superscripts {tt, ct, cc} correspond to the terms in the expression for AF^ (0) that 
involve the CKM factors {{VtlVtdf, (V;tKd)(^;6^td), (K^Kd)^} respectively. The subscript S 
denotes the contribution from the operator Qs, and the subscript 1/m denotes the terms 
that give the l/m^ corrections. The normalizing factor /C = {inml) / {3M^rr|BSo{xt)) and 
the value of (Am/F)^ may be taken from experiments: Xd = (Am/F)^ = 0.73 ± 0.03 0. 
The form of eq. (|58D can bring out important features of the dependence of AF^ on various 
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parameters, as we shall see below. This representation also has the advantage that within 
the leading term the CKM dependence cancels out and the value of Xd is available from 
experiments. 

A remark about the penguin contributions is in order. We only include the interference 
of the penguin operators C3...CQ with the leading operators Ci and C2. At the NLO, this 
approximation can be made consistent (in the sense of scheme independence) by counting 
the Wilson coefficients C^.-.C^ as of order as- These Wilson coefficients are modified at NLO 
through the mixing of Ci and C2 into C^-.-Cq. For Ci and C2 we use the complete NLO 
values. Since the contribution due to Cg starts only at the NLO level, we only have to use 
the LO value for that Wilson coefficient. We stress that if one uses the consistent NLO 
approximation just described, the corresponding LO approximation includes no penguin 
contributions and uses the LO values for Ci and C2. 

The choice of the 6-quark mass at LO is ambiguous (it may be taken to be the pole mass 
or the running mass at one or two loop level); we take it to be the running mass in the MS 
scheme to leading order in m?,. 

We use the following values of parameters to estimate Ar^: 

Mb^ = 5.28 GeV , = 4.8 GeV , m, = L4 GeV , 

mi,{mb) = 4.4 GeV , Tht{mi,) = 167 GeV . (59) 

To the NLO precision [we use here the NDR scheme to get r]B{mb) = 0.846 and include 



the NLO Wilson coefficients [14] and the corrections computed in eqs. (49), (51)], we get (in 
units of 10~^) 

/Ar\ Bs 1 
= n T.'^ -^ 4 Q.'^— 



, ^ , 0.73 + 4.93-^-1.38- 



cos/5 / Bs I { Bs 1 

^ ' 1.07 + 0.29-f + 0.02— + -0.005 - 0.021-f - 0.003— 



+^ (0.02 -0.06^ + ^ -0.01^ n . (60) 



Rt \ B B { B B 

»^(0.0.-0.06f.{-0.0l| 

Let us perform a conservative estimate of the error on the value of AF^^/Frf that we obtain 
here. The errors arise from the uncertainties in the values of the CKM parameters, the bag 
parameters and the mass of the b quark. There are also errors from the scale dependence, 
the breaking of the naive factorization approximation, and the neglected higher order terms 
in the z expansion. 

In the SM, we have 

cos = 1.03 ±0.08 , cos 2/?/i?^2 = 0.87 ±0.15 , (61) 



where we have taken the values of the CKM parameters from the global fit |T^. The leading 
term on the first line in ( |60l) is independent of the CKM elements. The quantity cos f3/Rt 
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is known to an accuracy of about 10% and appears in (pO) with a coefficient ~ 0.3 relative 
to the leading term. The quantity cos2/3/i?j, although known to only about 20%, appears 
with a very small coefficient (~ 10~^) as compared to the leading term in (|60|) . The net 
error due to the uncertainty in the CKM elements is thus approximately only 3%, i.e. about 
±0.1 X 10"^ 





LO 


A 


B 


C 


B 


1.0 


0.90 


0.83 


1.0 


Bs 


1.0 


0.75 


0.84 


1.0 


Ar,/r, 


6.3 X 10^^ 


1.9 X 10^^ 


2.6 X 10^^ 


2.8 X 10"^ 



Table 1: The numerical value of ATd/T for different values of the bag parameters. The 
column LO (C) shows the leading (next-to-leading) order result with factorization, i.e. B = 
Bs = 1. The values of the bag factors in column A are taken from and the ones in 
column B from the (preliminary) results in an unquenched {Nf = 2) lattice calculation by 
the JLQCD collaboration [g. 



We estimate the effect of the uncertainties in the bag factors by computing ( pOf ) with 
three sets of values of the bag parameters. The numerical results are as shown in Table |T[ 
From the table, and using the uncertainties on the values of the bag parameters as given in 
||16|| , we conservatively estimate the corresponding uncertainty in the value of AF^/r^ due to 
bag factors to be approximately ±0.5 x 10~^. The uncertainty in the value of rfif, = 4.4 ±0.2 
also leads to an error of ±0.5 x 10^'^. The uncertainty due to the scale fii dependence is 
estimated to be lii'^ x 10~^ (where fxi is varied between 2mf, and mh/2 following the common 
convention). The error due to the input value of Xd is 0.1 x 10~'^. 

The errors due to the breaking of the naive factorization assumption (which was made in 
the calculation of the matrix elements of the l/ruf, operators) are hard to quantify. Assuming 
an error of 30% in the R matrix elements (as in |T^), we estimate the error due to this source 
to be ±0.3 X 10-3. 

Table |I] also gives the LO value of ATd/Ta in the factorization approximation. We observe 
that the NLO corrections significantly decrease the value of ATd/Td as computed at LO, 
and that there effectively is no real (as/47r) suppression of the NLO contribution, as one 
naively expects. Therefore higher-order terms in the z expansion become important. While 
we estimate the error due the z expansion in the l/rrib and the penguin contributions from 
the terms in curly brackets in ( p(]| ) to be less than ±0.05 x 10^^, the issue of higher order 
terms in the NLO contribution ( |5B| ) is more subtle. We can write F21 in the form (see Sec. [4.1| 
for details) 

r2i(i?d) = -Ar[(\/;v;d)2[/(^,z)-2/(z,o)±/(o,o)] 

±2(v;;v;d)(v;*,\4d)[/(o,o) - fiz,o)] + (\/;,Vi,)V(o,o)l . (62) 
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Here in ( |55| ) we have included the complete NLO coefficient of (V^fel^d)^, which includes all 
the terms of the order in the hadronic matrix elements /. However, in order to calculate 
the corrections due to higher order terms in a complete NLO calculation is necessary. The 
contribution of these terms to [(/(O, 0) — j{z, z))l f{Q, 0)] can be computed to NLO precision 
using to be 

[(/(O, 0) - f{z, z))/fiO, 0)Ulo ^ 0.2 . (63) 

If we estimate the contribution to [(/(0,0) — f{z,0))/f{0,0)]NLo to be also of the same 
order, this results in the estimation of the net error in ATa/Td due to these terms to be 
±0.8 X 10-3. 

Our net estimation for the width difference is 

'^^) = (2.6t}J) X 10-=^ . (64) 

We have taken the central value to be the one obtained from the latest preliminary (un- 
quenched) results from lattice calculations [1^. The dominating theoretical errors are the 
scale dependence and the terms in F^-^'^ that correspond to the nontrivial z dependence 



of the function f{z,0) in (|6^) . To take care of the latter, a complete NLO calculation is 
definitely desirable. 

In the above calculations, we have used the expansion of in the form 



VlVr 



cb^cd 



ct 



ytbVtd, 



(65) 



^ -M{y:,Vuff{QM^ + 5fc 

where = 2[/(0, 0) - f{z, 0)]//(0, 0) and = [/(O, 0) - 2f{z, 0) + f{z, z)]l f{Q, 0). 

Following the suggestion in [|1310, we have also performed the expansion (and the error 
analysis) in the form 



F2i(S, 



-U{V:,Vtaff{z,z) 



ut 



^ub^ud 



^ -J^iV,lVtdffiz,z)il + 6fu) 



(66) 

where J^^t = 2[f{z, z) - f{z, 0)]/f{z, z) and T^u = [/(O, 0) - 2f{z, 0) + f{z, z)]/f{z, z). In 
this expansion the unknown NLO terms are suppressed by small CKM factors. This gives 
the width difference as 



AT, 



(s.olH) X 10- 



(67) 



where, as before, we use the latest preliminary (unquenched) results from lattice calculations 
for the bag parameters. The results of both (|6^) and (|67|) are consistent. The errors 



in (|67|) are smaller, but it should be noted that in both calculations the errors due to NLO 
terms are based on the assumption on the function /(z, 0) which we stated above. 
^We thank Uli Nierste for bringing this to our attention. 



14 



3 Measurement of AFd/Td 



It is not possible to find a final state to which the decay of i?^ involves only one of the decay 
widths r^, and Th- Indeed, since the Bd-Bd mixing phase (2/3) is large, the CP eigenstates 
are appreciably different from the lifetime eigenstates. The decay rate to a CP eigenstate 
therefore involves both the lifetimes. The semileptonic decays are flavor-tagging, and hence 
also involve both the lifetimes in equal proportion. 

We start by concentrating on the untagged measurements, i.e. the measurements in 
which the (Amt) oscillations are cancelled out. When the production asymmetry between 
Bd and Bd is zero (as is the case at the B factories), this corresponds to not having to 
determine whether the decaying meson was Bd or Bd- Restricting ourselves to untagged 
measurements is a way of getting rid of tagging inefficiencies and mistagging problems. At 
hadronic machines, a handle on the production asymmetry between Bd and Bd is necessary. 

In this section, we show explicitly that the time measurements of the decay of an untagged 
Bd to a single final state can only be sensitive to quadratic terms in /Wd/Td- This would im- 
ply that, for determining ^Td/Td using only one final state, the accuracy of the measurement 
needs to be (AP^^/r^)^ ~ 10"^. This indicates the necessity of combining measurements from 
two different final states to be sensitive to a quantity linear in ATd/Td ~ 0.3 x 10~^. We 
discuss three pairs of candidate channels for achieving this task. Finally, we point out the 
extent of systematic error in the conventional measurement of (3 due to the neglect of the 
width difference, and show how the tagged Bd J/ipKs mode can also measure AP^/P^ 
by itself. 

3.1 Quadratic sensitivity to ATd/Td of untagged measurements 

It is "common wisdom" that the time measurements, in general, are sensitive only quadrat- 
ically to AP^/Prf. Specific calculations (e.g. see 0]) also get results that can be clearly seen 
to obey this rule. Here, we give an explicit derivation of the general statement, pointing out 
the exact conditions under which the above statement is valid. Ways of getting around these 
conditions lead us to the decay modes that can provide measurements sensitive linearly to 



The non-oscillating part of the proper time distribution of the decay of Bd can be written 
in the most general form as 



The non-oscillating part can also be looked upon as the untagged measurement. 

For an isotropic decay, the only information available from the experiment is the time t. 



ATd/Td. 
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This information may be completely encoded in terms of the (infinitely many) time moments 

Jt-fit)dt 



(0 = 



Jf{t)dt 



(69) 



Expanding in powers of ATd/Td, we get 

in 



1 TT^ + C* 

2 V, 



(Ar,/r, 



(70) 



Defining the effective untagged lifetime as Tf, = (1 — available information 



(71) 



(p9|) is encoded in 

{t^) = n\inr [i + o [(Ar,/r,)2" 

Thus, when the accuracy of the lifetime measurement is less than (Ar^/r^)^, only the com- 
bination Tf, of Frf, AFrf and b may be measured through a single final state. This measurement 
is insensitive to b (to this order) and hence incapable of even discerning the presence of two 
distinct lifetimes (6 = and 6=1 would correspond to the presence of only a single lifetime 
involved in the decay.) In particular, in order to determine ATd/Td, the lifetime measure- 
ment through the semileptonic decay needs to be more accurate than {ATd/Tdy ~ 10~^. 
This task is beyond the capacity of the currently planned experiments. 

Combining time measurements from two different final states, however, can enable us to 
measure quantities linear in AFrf/r^. Indeed, for two final states with different values 6 (say 
6i and 62), we can measure 



_ ^ _^ 62_ - 61 ATd 



(Ar,/r,) 



(72) 



In the next subsections, we discuss pairs of decay channels that can measure this quantity 
([7^) that is linear in ATd/Td- 

3.2 Decay widths in semileptonic and CP-specific channels 

Let us first develop the formalism that will be applicable for all the decays that we shall 
consider below. When the width difference is taken into account, the decay rate of an initial 
Bd to a final state / is given as follows. Let Aj = {f\Bd), Af = {f\Bd), and 



A 



/ 



qAf_ 



(73) 



where p and q are as defined in (^). Using the CP- violating parameter 6'^ as defined in (|^ 
we get 

Q 



p 



(74) 



16 



The approximation here is vahd since we have l^"^! ~ \/S.Td/ ^md\ < 10~^. Henceforth, we 
shall only consider terms linear in 5°'. 

The decay rate of an initial tagged Bd or Bd to a final state / is given by f^: 

|2 



V{Bd{t)^f) =Uf\A^ 



1 + |A 



/I 



/I 



Ar t Ar t 

cosh — + cos(Amt) + A/^v sinh ^ + sin (Ami) 



1 + |A/ 



Ar t Ar t 

cosh ^ - A%1> cos(Amt) + ^af sinh — - ^^^^ sin(Amt) 



where the CP asymmetries are defined as 

|2 



A!: 



dir 
CP 



1- |A 



/I 



1 + |A 



|2 ' 



J mix 
•^CP 



2ImA 



/ 



/I 



1 + lA, 



and 



2ReA/ 

l + IA, 



|2 ' 



(75) 



(76) 



(77) 



and Mf is a time-independent normalization factor. 

In the case of semileptonic decays, / = {Di^u}, so that Af = and hence A/ = 0. The 
time evolution then becomes 

ATdt 



TiBd{t)-^f) (X e 



cosh ■ 



cos(Amt) 



cx e + e ^ ^* + oscillating terms 



-Tnt 



(78) 
(79) 



so that for semileptonic decays, we have bsL = 0. Note that 6 = is true for all self-tagging 
modes, so that all the arguments below for semileptonic modes hold true also for all the 
self-tagging decay modes. 

For the decays to CP eigenstates that proceed only through tree processes (and have zero 
or negligible penguin contribution), we have Aj = ±(1 — 6'^)e~'^^^ (the two signs and 
"— " correspond to CP-even and CP-odd final states respectively). Then ([75|) gives 



r{Bd{t)^f) (X e-^^' 
oc e"^^* 



cosh — r— =F cos(2/9) sinh — ± sin(2/3) sin (Amt) 



2 ' ' ' 2 

l±cos(2/5))+e"^^*(lTcos(2/3))+ oscillating terms , (81 



0) 



where we have neglected the small corrections due to S'^. Thus, for CP eigenstates, we have 
bcp+ = +cos(2/5) and bcp- = — cos(2/?). 

The ratio between the two lifetimes tcp± and tsl is then 



TCP± 



(82) 



The measurement of these two lifetimes should be able to give us a value of lAr^l, since 
I cos(2/3)| will already be known to a good accuracy by that time. 
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Note that it is also possible to measure the ratio of the lifetimes tcp- and tcp+- 

l£^ = l + cosm^ + oUAT,/T,f] . (83) 

Although the deviation of the ratio from 1.0 in this case is larger by a factor of 2, using the 
effective semileptonic lifetime instead of one of the CP eigenstates would still be the favoured 
method. This is because the CP specific decay modes of Bd (e.g. J/ipKs(L), D^D^) have 
smaller branching ratios than the semileptonic modes. In addition, the "semileptonic" data 
sample may be enhanced by including the self-tagging decay modes (e.g. that 
also have large branching ratios. After 5 years of LHC, we should have about 5 x 10^ events 
of J/ipKsi whereas the number of semileptonic decays at LHCb alone that will be directly 
useful in the lifetime measurements is expected to be more than 10^ per year, even with 
conservative estimates of efficiencies. 

3.3 Transversity angle distribution in J /ipK* 

The decays B^, VV (where VV is a flavour-blind flnal state consisting of two vector 
mesons) take place both through CP-even and CP-odd channels. Since the angular informa- 
tion is available here in addition to the time information, these decay modes are not subject 



to the constraints of the theorem in Sec. |3.1| , and quantities sensitive linearly to Ar^/Ff^ can 
be obtained through a single final state. This cancels out many systematic uncertainties, 
and hence these modes can be extremely useful as long as the direct CP violation is negli- 
gible, and we can disentangle the CP-even and CP-odd final states from each other. This 
separation can indeed be achieved through the transversity angle distribution ([|T^-[|1^]). 



We illustrate the procedure with the example of B^ J/tlj{i~^i~)K*{Ks'n'^). The most 
general amplitude for the decay B J/ipK* is given in terms of the polarizations ej/^, ex* 
of the two vector mesons: 

A{Bd - J/i^K*) = A, (1^) e%e*^. - ^ e}J^ ■ ef. - e}/^ x e^,. ■ p , (84) 

where Ek* is the energy of the K* and p the unit vector in the direction of K* in the J/ip 
rest frame. The superscripts L and T represent the longitudinal and transverse components 
respectively. Since the direct CP violation in this mode is negligible, the amplitudes Aq and 
A\\ are CP-even, whereas A_[_ is CP-odd. Let us define the angles as follows. Let the x axis 
be the direction of K* in the J/ip rest frame, and the z axis be perpendicular to the decay 
plane of K* — > Kstt^, with the positive y direction chosen such that py{Ks) > 0. Then we 
define {6, (p) as the decay direction of in the J/ip rest frame and ip as the angle made by 
Ks with the x axis in the K* rest frame. 

Here 6 is the transversity angle, i.e. the angular distribution in 6 can separate CP-even 
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and CP-odd components of the final state. The angular distribution is given by 



dcosO 



+ cos'^) + ^\A4t)\\m^e (85) 



where = |^o(i)P + is the CP-even component and = \A±{t)\^ the 

CP-odd one. These two components can be separated from the angular distribution ( p5| ) 
through a likelihood fit or through the method of angular moments pT|f|. 
The time evolutions of the CP-even and CP-odd components are given by 



\A-it)\' 



A+(0)|' cos'p e-^^' + sin' f3 e~' "' + e-' "' sm{AMdt)sm{2(3) 



|A_(0)|" sin' /5 e"^^* + cos' /3 e 



-Fnt 



e"^^* sin(AMdt) sin(2/?) 



(86) 
(87) 



These are the same as the time evolutions in (0). The difference in the untagged lifetimes 
of the two components, 



TCP- 

TCP+ 



AT 

i + cos(2/?)— ^ + o (Ar,/r,) 



(88) 



is linear in the lifetime difference AF^. 

The disentanglement of the CP-even and CP-odd components from the angular distribu- 



tion is a statistically efficient process ||2T|. In fact, in the Bg system, the angular distribution 
of Bg J /ip{£'^ £~)(f){K^ K^) can be used for determining the lifetime difference AF^, and 
is the preferred mode for measuring this quantity. 

The mode J/ipK* suffers from the presence of a vr^ in the final state, which may be missed 
by the detector, thus introducing a source of systematic error that needs to be minimized. 



3.4 Untagged asymmetry between B J/ipKs and J/ipKi 

Two of the decay modes of Bd that have been well explored experimentally (because of 
their usefulness in measuring (3) are B J/ipKs and J/ipKi- Here we show that the 
time-dependent asymmetry between the decay rates of these modes is a quantity linear in 
ATd/Tii, and therefore within the domain of experimental feasibility. 
Let us define 

A{Bd ^ J/^Ks) = As, A{Bd ^ J/^Ks) = As, 
A{Bd ^ J/^jKl) = Al, A{Bd ^ J/^Kl) = Al, 

so that using 

\Ks) = {l + e)\K') + {l-e )\K') , \K = {I + e)\K') - {I - e)\K') , (89) 

^In we suggested to use the CP-odd-CP-even interference in the decay B — > J /ipK* to measure the 
value of APd/Pd. However, it involves tagged measurements in addition to two- or three-angle distributions, 
and hence is not as attractive as the untagged measurements described here. 
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we can write (with the phase convention Arg(g/p) = 0) 



As =Al= Ae'^(l + ape^V^^)(l + e) , 

As =-Al= Ae-'f'{l + ape'e-'^^){l-e) , (90) 

where ape^^e*^"^ is the ratio of contributions that involve the CKM factors V*^Vcs and V^lVtd 
respectively. The latter contribution (penguin) is highly suppressed with respect to the 
former one (tree): the value of is less than a percent. Here 6 is the strong phase and 
A7 = ATg{VtlVts/V*,,Vcd) ^ -0.015 in the SM. From (|7|), (]7|) and (^, we get 

\s = -\L^-il-5'^)e~^'^il-2e-2ismA-f ape'^) = -e-^'^{l-2e) , (91) 

where e is an effective complex parameter that absorbs all the small theoretical uncertainties. 

When the production asymmetry between and 13^ is zero (as is the case at the B 
factories), the untagged rate of decay is 



(92) 



T[B^n^J/i^Ks{KL)] ^ Af\As\\l-2Re{e))e-^'^' X 

cosh + ^Ar, sinh 

The only difference between the decay to Ks and that to Kl is the sign of ^Ar: 

AAr{Ks) = -AMKl) = cos(2/3) - 2 Im(e) sin(2/3) . (93) 

The untagged time-dependent asymmetry between Bun — > J/i^Ks and Kl is 

.fr. r.. rjBunit) ^ J/^jjKs) - r(i?„„(t) ^ J/^Kl) 

^ T{Bun{t)^J/^Ks)+T{Bun{t)^J/^KL) ^ ' 

= cos(2/3) tanh [1-2 Im(e) tan(2/3)] (95) 

^ cos(2/3) tanh (^^^) . (96) 

Thus, the measurement of this asymmetry will enable us to determine lAF^I, given sufficient 
statistics and a measurement of sin 2(3. 

The factor limiting the accuracy of the above asymmetry is the measurement of r(i?u„(t) — 
J/iIjKi). At the B factories, Kl may be detected through its hadronic interactions in the 
calorimeter, and though its energy is poorly measured, the corresponding detection of J/ip 
can help in reducing the background. However, the number of events available may be too 
small for an accurate measurement. In the hadronic machines this decay has a high back- 
ground, so the systematic errors in the measurement may be too large for this method to be 
of practical use. 
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3.5 Effect on the measurement of sm(2/3) 



The time- dependent CP asymmetry measured through the "gold-plated" mode J / ipKs 



IS 



21 H 



A, 



CP 



T[Bdit)^ J/^PKs]+T[Bdit) 
sin(Amrft) sin(2/3) , 



(97) 
(98) 



which is vahd when the lifetime difference, the direct CP violation, and the mixing in the 
neutral K mesons is neglected. As the accuracy of this measurement increases, the correc- 
tions due to these factors will need to be taken into account. Keeping only linear terms in 
the small quantities e and AP^, we get 



Acp = sin(Amt) sin(2/3) 



1 - sinh (^^Y~ ) 



+2Re{e) [-1 + sin^(2/5) sin^(Amt) - cos(Amt) 
+21m{e) cos(2^) sm{Amt) . 



(99) 

(100) 
(101) 



The first term in (^9]) represents the standard approximation used (|98|) and the correction due 
to the lifetime difference AP^. The rest of the terms [( |100| ) and (|101|) ] include corrections due 
to the CP violation in B B and K-K mixings, which are of the same order as AP^/ Gammad- 
In the future experiments that aim to measure /? to an accuracy of 0.005 the correction 
terms need to be taken into account. The corrections due to e and AP^ will form a major part 
of the systematic error, which can be taken care of by a simultaneous fit to sin(2/5), AP^ and 
e. The BaBar collaboration gives the bound on the coefficient of cos(Amt) in (|100|) , while 



neglecting the other correction terms [Q. When the measurements are accurate enough to 
measure the cos(Amt) term, the rest of the terms would also have come within the domain 
of measurability. For a correct treatment at this level of accuracy, the complete expression 
for Acp above (P ^101|) needs to be used. 



3.6 Tagged measurements 

Until now, we have discussed only the untagged measurements. Taking into account the 
oscillating part of the time evolution of the decay rate, we have the decay rate in general as 



g[t) = f{t) + Ce-'^'sm{Amt + $) 



(102) 



where f{t) is the untagged decay rate as defined in (pSj) , C a constant and $ a phase. The 
lifetime of the oscillating part is an additional lifetime measurement, which opens up the 
possibihty of being able to determine APd/P^ through only one final state (and without 
angular distributions as in Sec. p.3|) . 
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In the case of the semileptonic decays, this strategy fails since the semileptonic width 
measured with the untagged sample is 

so that 

TsL/Td = l + 0[{ATjTdy] . (104) 

Thus the semileptonic decays would provide sensitivity only to quadratic terms in ATd/Td, 
even if it were possible to use the tagged measurements efficiently. 

However, the untagged lifetime measured through the decay to a CP eigenstate is 

rcP± « (^1 T — , (105) 

SO that it differs from the lifetime of the oscillating part (tc; = l/F^) by terms linear in 
ATfi/Ta. Thus, the tagged measurements of a CP-even or CP-odd final state {D^D^, 
J/ipKs, J/ipK^, etc.) can measure ATa/Td by themselves. 

The mistag fraction is the main limiting factor on the accuracy of this measurement, and 
the tagging efficiency limits the number of events available. It is indeed possible that the 
measurement through the semileptonic decays will be more accurate than that through the 
oscillating part of the CP-specific final state. This then reduces to the method suggested in 



Sec. |3.2| . For further experimental details on a tagged measurement of ATd/Td we refer the 



reader to reference 12^ 



4 Lifetime differences in Bg and Bd systems 

The calculations of the lifetime difference in Bd (as performed here) and in the Bs system 
(as in [^, ^) run along similar lines. However, there are some subtle differences involved, due 
to the values of the different CKM elements involved, which have significant consequences. 
In particular, whereas the upper bound on the value of AP^ (including the effects of new 
physics) is the value of Ars(SM) ||26[, such an upper bound on AP^ can be established only 
under certain conditions and involves a multiplicative factor in addition to AP^(SM). Also, 
whereas the difference in lifetimes of CP-specific final states in the Bg system cannot resolve 
the discrete ambiguity in the Bs-Bg mixing phase, the corresponding measurement in the 
Bd system can resolve the discrete ambiguity in the Bd-Bd mixing phase. Let us elaborate 
on these two differences in this section. 
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4.1 Upper bounds on Ar(^(s) in the presence of new physics 

For convenience, let us define 6^ = Arg(r2i)g, $q = Arg(M2i)5, where q G {d, s}. Then we 
can write 

Ar, = -2|r2i|gcos(e,-$g) . (loe) 

Since the contribution to comes only from tree diagrams, we expect the effect of new 
physics on this quantity to be very small. We therefore take |r2i|5 and Qg to be unaffected 
by new physics. On the other hand, the mixing phase $q appears from loop diagrams and 
can therefore be very sensitive to new physics. 

Let us first consider the Bs system. Here r2i may be written in the form 

121(5.) = -Af[iV^V,s)'fiz,z) + 2iV:,V,,)iV:,Vus)fiz,0) + iV:,Vus)'fiO,0)] (107) 

where A/" is a positive normalization constant and f{x,y) are the hadronic factors that do 
not depend on the CKM matrix elements. In the limit z = ml/ ml — > 0, we get f{z,z) = 
f{z,0) = /(0,0). Even with z ^ 0, we expect all the /'s to have similar magnitudes. 
On the other hand, the CKM elements involved in (|107|) obey the hierarchy (V^^Ks)^ ~ 
A' , {V:,Vcs){V:,Vus) ~ , {V:,VusY ~ A^. The term involving (V;*,^,)^ then dominates in 
(|107|), and we can write 



121(5.) = -AfiV:,V,s)'fiz,z)[l + OiX')] . (108) 
Since the /'s are real positive functions, we have 6. ~ vr + Arg(V^'^\4.)^. Then, 

Ar. = 2|r2i|.cos[Arg(v;;Ks)'-$.] • (io9) 

In SM, $s = Arg(l^^Vts)^, therefore the argument of the cosine term in (|109|) is given by 
Arg[{V;^V^sf/{VtlVtsf] = -2A7 ^ 0.03. Thus in SM, we have 

Ar,(SM) = 2|r2i|sCos(2A7) . (110) 

The effect of new physics on AF^ can then be bounded by giving an upper bound on AF^: 

AF, < ^ AF.(SM) . (Ill) 

cos(2A7) 



Thus, the value of AF^ can only decrease in the presence of new physics 

In the case of the Bd system, the situation is slightly different. As in the Bg case, we can 
write 

r2i(5,) = -^^[{v:M'f{z,z) + 2{v:Miv:,Vud)fiz,o) + (kT^k^) V(o, o)] (112) 

where the normalizing factor Af and the hadronic factors / are the same as in the Bg case 
in the limit of the U-spin symmetry (see [27|) and therefore have similar magnitudes. The 
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CKM elements involved in ( |112D do not obey a hierarchy similar to the Bg case: instead 
we have ~ (KfcKd)!^:^^^) ~ (KT^Kd)' ~ Then no single term in can 

dominate. We can, however, use the unitarity of the CKM matrix^ to rearrange ( |112|) in the 
form 

T2iiB,) = -^^[(v:,v,,Y[f{z,z)-2f{z,o) + fio,o)] 

+2iV:,V,,)iV*Vt,)[fiO,0) - /(z,0)] + iV*Vt,)'fiO,0)] . (113) 

Note that in the limit of z — 0, all the factors / are identical and hence the coefficients 
of (y*fycdY and (y*^Vcd)(y^lVtd) vanish. The last term in (|TT^) is then left over as the 
dominating one, and we get 

r2i(5d) ^ -^fiv,lVtd)'fio, 0) , (ii4) 

thus giving ~ tt + Arg(V^^Vjd)^. However, the finite value of z ~ 0.1 may give large 
corrections to this value. 

We have already shown that ( p. 14]) may be written in the forms (|65D or (|66D. From this, 
we get 6rf = TT + Aigiy^lVidY + Arg(l + 6fi), where i G {u, c}. Using (|106|) , we then have 

ATd ^ 2\T2i\dCos[ATg{V;,Vuf -^d + Arg(l + 6f,)] . (115) 

In SM, $rf = AigiV^lVtdf, so that 

Arrf(SM)^2|r2iUcos[Arg(l + 5/,)] , (116) 

and an upper bound for AF^ can be written in terms of AFd(SM) as 

cos[Arg(l + 6fi}\ 

We can calculate the bound (|117|) in terms of the extent of the higher order NLO corrections. 



Estimating \J-'ct\ < 0.4 (corresponding to the error analysis in Sec. |2.4|) , we get |Arg(l + 
Sfu)\ < 0.6, so that we have the bound ATd < 1-2 ATd{SM). Note that this bound is valid 
only in the range of J-'ct estimated above. A complete NLO calculation will be able to give 
a stronger bound. 

Thus in the case of the Bd system, we have an upper bound analogous to the one in the 
Bg system only under certain conditions. Moreover, the reasons behind the existence of these 



^ We note that this assumption of the unitarity for a three-generation CKM matrix is quite general, 
because most popular new physics models, including supersymmetric models, preserve the three-generation 
CKM unitarity. The present CKM values, constrained from various experiments, are completely consistent 
with the unitarity for the three-generation CKM matrix. Moreover, one can show that the non-unitary effects 
within the three-generation CKM, which can stem from the fourth generation or E(6)-inspired models with 
one singlet down- type quark, are < A^, once we assume a Wolfenstein-type hierarchical structure for the 
extended CKM matrix. 
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two upper bounds differ. Whereas in the Bg case it follows directly from the hierarchy in 
the CKM elements, in the case it depends on the values of the hadronic terms. Note that 
whereas unitarity was not needed in the Bs case, the assumption that (r2i)5 is unaffected 
by new physics is required in both the cases. 



4.2 Mixing phase: new physics and discrete ambiguity 

The Bd-Bd mixing phase is efficiently measured through the decay modes J/ipKs and 
J/ipKi. If we take the new physics effects into account, the time-dependent asymmetry 
(13) is Acp = -sin(AMdt)sin($d), which reduces to (H) in the SM, where $rf = -2(3. 
The measurement of sin($(i) still allows for a discrete ambiguity $rf ^ vr — Whenever 
a discrete ambiguity in P is referred to (/5 ^ n/2 — P) in this paper (or in the literature), 
strictly speaking we are talking about the discrete ambiguity $d ^ — ^d- In this section, 
we shall use the notation $rf instead of 2(3 in order to illustrate the comparison with the 
corresponding quantities in the Bg system. 

Getting rid of the above discrete ambiguity is a way of uncovering a possible signal of 
new physic^. Ways to get rid of this ambiguity have been suggested in literature, using 
the comparison of CP asymmetries in J/ipKs and vrvr |2^, time dependent CP asymmetries 



in Bg — > pKs and in Bg — > ttK, KK |^T|, angular distributions and U-spin symmetry 
arguments [0, or cascade decays B ^ D ^ K |Q. The measurement of $d through the 



measurements involving AF^^ is unique in the sense that it uses only untagged measurements. 

In Sections |3.2| and |3.3| , we have seen that the ratio of two effective lifetimes can enable 
us to measure the quantity /Wohs{d) = cos{2P)ATd/T. In the presence of new physics, this 



quantity is in fact (see eq. |106|) 



In SM, we get 



ATobsid) = -2i\T2i\d/Td) cos(<l>d) cos(0d - $rf) . (118) 



Aro,,(,)(SM) = 2(|r2i|d/r,) cos(2/3) cos[Arg(l + 5fi)] . (119) 



If \Sfi\ < 1.0, we have cos[Arg(l + 6f)] > (in fact, from the fit in |T5[ and our error 
estimates, we have cos[Arg(l + 5/„)] > 0.8). Then Aro;,s(d)(SM) is predicted to be positive. 
New physics is not expected to affect 6^, but it may affect $d in such a way as to make the 
combination cos($ri) cos(0rf — $d) change sign. A negative sign of ATobs{d) is a clear signal 
of such new physics. 

^ In SM, the value of must match with the phase of the b d penguin. However, the direct 
measurement of the latter phase is not theoretically clean so the preferred way is to compare the 
measured value of with the value of 2/3 determined through an "unitarity fit" for all the CKM parameters 
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In addition, if can be determined independently of the mixing in the Bd systemQ 
then measuring ATobs{d) and using ( |118| ) gives us two solutions for in general: $1^ and 
$2d such that tan($i(i + $2d) = tan(6d). As long as tan(6(i) 7^ (as the unitarity fit 
suggests), at most one of the solutions will correspond to the value of sm{^ a) j/ipKs obtained 
through AcpiJ/i'Ks). If none of $^1 or $^2 corresponds to the value of sm{^d)j/^pKsy there 
definitely is new physics in mixing If one of $1^ or $2^ matches with sin($d) j/V'ii's) 
it gives the actual value of and thus resolves the discrete ambiguity in principle. In 
practice, this requires the knowledge of r2i(d) theoretically to a high precision and having 
to measure ATobs(d) to sufficient accuracy to be able to distinguish between $1^ and $2d- A 
complete NLO calculation is needed for the former. The latter may be achieved at the LHC 
using the effective lifetimes of decays to semileptonic final states and to J/ipKs- 

Let us contrast this case with that in the Bg system and show these features above are 
unique to the B^ system: The corresponding time-dependent asymmetry in the Bg system is 
measured through the modes J/4'4' o^' J/^V^'\ which give the value of sin($<j), and therefore 
leave the discrete ambiguity $5 ^ vr — $s unresolved. The ratio of two effective lifetimes in 
the Bg system can enable us to measure the quantity 

ATobsis) = cos($,)Ar,/r 

= -2|r2i|,/r,cos(<i>,)cos(e,-<i>,) . (120) 

Since Qs ^ 71 + Arg(l/^Vcs)^ ~ vr, we have 

Ar,bs{s) ^ 2|r2i|./r,cos2($,) . (121) 

This measurement thus still has the same discrete ambiguity $s ^ — $s as in the J/ip(j) 
or case, and the discrete ambiguity in the Bg system is not resolved. 

5 Summary and conclusions 

It has been known for many years that the B^ system is a particularly good place to test the 
standard model explanation of CP violation through the unitary CKM matrix. The phase 
2(3 involved in the Bd — Bd mixing is large, and hence the CP violation is expected to be 
larger in the Bd system in general, as compared to the K or the Bg system. This feature 
has already been exploited in various methods for extracting a, (3 and 7, the angles of the 
unitarity triangle, by measuring CP-violating rate asymmetries in the decays of neutral Bd 
mesons to a variety of final states. In particular, the precise measurement of sin (2/5) from 

""^"For example, if we assume that new physics does not affect the mixing in the neutral K system, the fit 
to IKib/Vcbl and ck (a part of the unitarity fit) determines Arg(V'j^Vt(i) independent of any mixing in the B 
system. This gives us a measurement of 0^ = tt + Piig[V^fytdf + Arg{l + Sfi) where Arg{l + Sfi) will be 
known to a good accuracy with a complete NLO calculation. 

Or any assumption that goes into the determination of Qd is in question. 
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the theoretically clean decay modes Bd{t) J /i]jKs{Kl) is a test of the SM, as well as the 
opportunity to search for the presence of physics beyond the standard model. 

The two mass eigenstates of the neutral system — Eh and El — have slightly 
different lifetimes: the lifetime difference is less than a percent. At the present accuracy of 
measurements, this lifetime difference AF^^ can well be ignored. As a result, the measurement 
and the phenomenology of AF^^ has been neglected so far, as compared to the lifetime 
difference in the Es system for example. However, with the possibility of experiments with 
high time resolution and high statistics, such as the electronic asymmetric E factories of 
BaBar, BELLE, and hadronic E factories of CDF, LHC and BTeV, this quantity starts 
becoming more and more relevant. 

Taking the effect of AFj^ into account is important in two aspects. On one hand, it affects 
the accurate measurements of crucial quantities like the CKM phase (5 and therefore must 
be measured in order to estimate and correct the error due to it. On the other hand, the 
measurement of AF^ can lead to clear signal for new physics. 

Thus in addition to being the measurement of a well-defined physical quantity which can 
be compared with the theoretical prediction, the value of AF^ is important for getting a firm 
grip on our understanding of CP violation. It is therefore worthwhile to have a look at this 
quantity and make a realistic estimation of the possibility of its measurement, as we do in 
this paper. 

Wc estimate AF^/F^ including l/mj, contributions and part of the next-to-leading or- 
der QCD corrections. We find that adding the latter corrections decreases the value of 
AFd/F(i computed at the leading order by almost a factor of two. We get the final result 
as AF(i/Fd = (2.61|;q) x 10"'^, or using another expansion of the NLO QCD corrections 
AFrf/F(i = (3.0l?;4) x 10^^, where for the central value we have used the preliminary values 
for the bag factors from the JLQCD collaboration. In the error estimation, wc take into 
account the errors from the uncertainties in the values of the CKM parameters, the bag 
parameters, the mass of the b quark, and the measured value of Xd- The major sources of 
error are the scale dependence, the breaking of the quark-hadron duality, and the error due 
to missing terms in the NLO contribution. We show that a complete NLO calculation is 
desirable. 

The most obvious way of trying to measure the lifetime difference is through the semilep- 
tonic decays, however it runs into major difficulties. If only the non-oscillating (untagged) 
part of the time evolution of the decay is considered, we indeed have a combination of two 
exponential decays with different lifetimes. However, as we show in this paper, there is no 
observable quantity here that is linear in AF^i/Fd. The time measurements allow us to de- 
termine the quantity tsl = d)[^ -\- (^{^^ dl^ df\- This decay mode is thus sensitive only 
to quantities quandratic in AF^/F^. So this method would involve measuring a quantity as 
small as (AF^/Fd)^ ~ 10~^, which is not practical. The lifetime of the oscillating part is 
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also l/Frf, so adding the information from the oscillating part of the time evolution does not 
help at all. This problem arises for all self-tagging decays. Therefore, though self-tagging 
decays of have significant branching ratios, they cannot by themselves be expected to 
give a measurement of Ar^/F^. 

The time evolutions of Bd decaying into CP eigenstates also involve both the lifetimes, 
since the Bd — B^ mixing phase (2/9) is large, which implies that the CP eigenstates are 
appreciably different from the lifetime eigenstates. As a result, it is not possible to find a 
final state to which the decay of B^, involves only one of the decay widths and Tu. The 
non-oscillating part of the time evolution of decays to CP eigenstates gives a quantity rcp± = 
(l/rd)[l ± {cos{2(5)/2)ATd/Td + 0{M:d/Td)% but the quantities and AF^ cannot be 
separately determined through this measurement, and sensitivity to (AF^/Fd)^ is necessary. 
(Indeed, we explicitly prove a general theorem that shows that, for isotropic decays of Bd to 
any final state, the untagged measurements can only be sensitive to (AFji/F^)^.) 

The oscillating part of the time evolution to CP eigenstates has a hfetime l/F^ (to 
an accuracy of 0(AFd/Fd)^). Therefore, if this hfetime is measured accurately, it can be 
combined with the measurement of tcp± through the untagged part to get a measurement 
linear in AFd/F^. However, the need for tagging, and consequent mistagging errors, reduce 
the efficiency of this method. 

A viable option, perhaps the most efficient among the ones considered here, is to compare 
the measurements of the untagged lifetimes tsl and tcp±. Since tsl is in fact the lifetime 
for all self-tagging decays, and the branching ratios for self-tagging decays of Bd are much 
larger than the decays to CP eigenstates, we expect that the most useful combination will be 
the measurement of tsl through self-tagging decays and that of tcpj^ through Bd J /ipKs- 

The untagged asymmetry between Bd J/ipKs and Bd — ^ J/i/jKl is a particular case 
of using the combination of measurements of tcp^ and tcp_, which we analyze in detail. The 
effects of CP violation in the mixing and decay of Bd, as well as the indirect CP violation 
in the K system has been taken into account. The usefulness of this method is limited by 
the uncertainties in the measurement of r{Bd — > J/ipKi). 

Since the theorem referred to above — about a single untagged decay being sensitive 
only to (AFrf/Frf)^ — applies only to isotropic decays, decays of the type B — > VV can still 
be used by themselves to determine quantities linear in AFd/F^. A promising example is 
Bd t't') K*{-^ KsTT^). The CP-odd and CP-even components in the final state 

can be disentangled through the transversity angle distribution, and both tcp+ and tcp_ can 
be determined through the same decay. Since there is only one final state, many systematic 
errors are reduced. The only undesirable feature of this decay mode is the presence of 7r° 
in the final state, which may be missed, especially in the hadronic machines. The three 
angle distribution of the same decay mode can also be used to obtain AF^/F^ through the 
interference between CP-even and CP-odd final states. The three angle method is however 
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not as efficient as tfie single angle distribution, since one has to use tagged decays and more 
number of parameters need to be fitted. 

We also point out the interlinked nature of the accurate measurements of f3 and Ar^/rd 
through the conventional gold-plated decay. In the future experiments that aim to measure 
(3 to an accuracy of 0.005 or better, the corrections due to AT^ will form the major part of 
the systematic error, which can be taken care of by a simultaneous fit to sm{2 j3), AT ^ and 
an effective parameter e that comes from a combination of CP violation in mixing in the B^, 
and K system. 

All the combinations of untagged decay modes discussed here involve measuring the 
quantity ATobs{d) oc (cos(2/9)/2)Ard/r(i, wherein the value of AFd/Fd also depends on p. 
The complete dependence on (3 is of the form cos(2/3) cos(0^+2/3) . The sign of this quantity is 
positive in SM, but may be changed in the presence of new physics. There is thus a potential 
for detecting physics beyond the SM. Moreover, if the value of is known independently 
of B mixing, then since Arobs{d) is not invariant under /3 <-> 7r/2 — /3, the discrete ambiguity 
in P is resolved in principle. Note that this feature is unique to the B^ system — in the Bg 
system for example, Aro6s(s) does not help in resolving the corresponding discrete ambiguity 
in the Bg-Bg mixing phase. 

It is known that, if (r2i)s is unaffected by new physics, then the value of AF^ in the Bg 
system is bounded from above by its value as calculated in the SM. In the B^ system, this 
statement does not strictly hold true. However, if (r2i)ci is unaffected by new physics and 
the unitarity of the 3x3 CKM matrix holds, then an upper bound on the value of AF^^ may 
be found. In the absence of a complete NLO calculation, however, this bound is a weak one. 

With the high statistics and accurate time resolution of the upcoming experiments, the 
measurement of AF^ seems to be in the domain of measurability. And given the rich phe- 
nomenology that comes with it, it is certainly a worthwhile endeavor. 
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